Common Origin of Quantum Regression 
and Quantum Fluctuation Dissipation Theorems 

P. Shiktorov, E. Starikov, V. Gruzinskis 
Semiconductor Physics Institute, A. Gostauto 11, 2600 Vilnius, Lithuania 

L. Reggiani 

Dipartimento di Ingegneria dell' Innovazione, 
Istituto Nazionale di Fisica della Materia, 
Universita di Lecce, Via Arnesano s/n, 73100 Lecce, Italy 

L. Varani and J.C. Vaissiere 

Centre d'Electronique et de Micro-optoelectronique de Montpellier 
(CNRS UMR 5507), Universite Montpellier II, 34095 Montpellier Cedex 5, France 

Abstract. 

It is shown that the quantum fluctuation dissipation theorem can be considered as a 
mathematical formulation in the spectral representation of Onsager hypothesis on the 
regression of fluctuations in physical systems. It is shown that the quantum fluctuation 
dissipation theorem can be generalized to an arbitrary stationary state. 
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Introduction. Under thermal equilibrium conditions the behavior of fluctuations of 
macroscopic observables of a physical system is governed by relationships which are 
formulated usually in terms of the regression theorem (the so called Onsager hypothesis 
[1]) and the fluctuation-dissipation theorem [2-4]. The former pertains to the time 
domain and states that the relaxation of a correlation of fluctuations is described by the 
same law governing the irreversible processes of the observable quantity itself. The latter 
pertains to the frequency domain and interrelates in some universal way the spectral 
characteristics of fluctuations and linear response (i.e. dissipation) of an observable 
of the physical system. In the framework of a classical approach, the two theorems 
complement each other providing a closed description of fluctuations and linear response 
in thermal equilibrium. By contrast, within the more general quantum approach there 
appears a conflict between these two theorems. This conflict is usually interpreted 
as a violation of the quantum regression theorem (QRT) (see, for example, Refs. [5- 
7]). In the most evident form such a violation is demonstrated in Ref. [7], where the 
conclusion "there is no quantum regression theorem" is announced. The proof of the 
general character of such a statement is based on the fact that the violation of QRT 
follows from the quantum fluctuation dissipation theorem (QFDT). A proof that QRT is 
valid independently of QFDT was given by Lax [8] on the basis of the general principles 
of quantum statistics (see also Refs. [9-11]). Since from Ref. [7] it is claimed that 
QFDT contradicts the validity of QRT, we argue that the origin of such a conflict is 
related with QFDT and its interpretation. 

The aim of this letter is to establish the origin of such a conflict and to provide 
arguments to reconcile the two theorems. For this sake, the conditions necessary for 
QRT to be fullfilled are firstly considered in the framework of the general formalism of 
linear response functions. Then, the equivalence between QFDT and QRT is proven 
by revisiting the derivation of QFDT. Finally, an extension of QFDT to an arbitrary 
stationary state which can be far from that of thermal equilibrium is presented. 

The Onsager regression theorem. Let us define formally a law which will describe 
the time evolution of perturbations of some observable physical quantity x(t). It is 
assumed that the response of the observable, induced by some external perturbation 
described by a force /(£), is linear and can be represented as a convolution integral: 



Here a x (r) is a real function of time determined by the state and the properties of the 
physical system under test. The response described by Eq. (1) satisfies the causality 
principle, i.e. a x (r) = at r < 0. Therefore, one always can introduce a linear 
integro-differential operator A t {...}, such that x(t) given by Eq. (1) is a solution of the 
equation: 



Due to the linearity of A t {...}, its eigenfunctions are the harmonic functions exp(-iujt) , 
and the corresponding eigenvalues A(u) are related with the generalized susceptibility 
of the physical system a x (u) by: 
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(1) 



A t {x(t)} = f(t) 



(2) 




(3) 
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with a x (r) being the correspondent response function. By using the reciprocal Fourier 
transformation, from Eq. (3) it is easy to show that a x (r) satisfies the equation: 

A t {a x (r)} = S(t) (4) 

i.e., it describes the response of x(t) at r > induced by an impulsive perturbation at 
r = 0. 

The mathematical formulation of the Onsager hypothesis, i.e. the so called regres- 
sion theorem is that the correlation function of fluctuations of the observable x, C xx (t), 
must satisfy at r > the same equation as a x (r) at r > 0, i.e. 

A t {C xx (r)} = for r>0 (5) 

As a consequence of this theorem, the time dependence of both a x (r) and C xx (t) is 
described by the same set of eigen-solutions exp(—iuJut) obtained from the homogeneous 
equation: A t {exp(— iu>nt)} = 0. The resonant frequencies ur satisfy the dispersion 
equation A(ujr) = 0, i.e. in the complex frequency plane u = uj' + ioj" they are the 
poles of the susceptibility. Due to the causality principle, all poles of a x (u) are placed 
in the lower half-plane [3], i.e. uj r < 0, what provides the relaxation process. 

The response function formalism allows us to describe macroscopic properties of 
the physical system by considering merely the analytical properties of a x (u>). Let us 
consider the reasons leading to the violation of QRT by using the usual interrelation 
between the spectral density of the fluctuating observable, S xx (u>), and that of the force, 
Sff(w) [3], 

S xx {uj) = a x {ijj)a x {-Lo)Sff{uj) (6) 

In the framework of linear response formalism, the correlation function of observable 
fluctuations obeys the following integro-differential equation: 

1 f°° 

A t {C xx (r)} = — / a x (-u)Sf f (u)exp(-iUT)du (7) 

^ J-oo 

As follows from Eqs. (6) and (7), the fluctuations are governed by the properties of 
internal and external interactions. The former is associated with the resonances of 
the susceptibility a x (u>) and the latter with the frequency dependence of Sff(u>). As 
follows from Eq. (7), the internal interactions can not be a source of violation of the 
regression hypothesis since, due to the causality principle, a x (uj) has no poles in the 
upper half-plane, and, hence, a x (—u>), has no poles in the lower half-plane. Therefore, 
the right hand side of Eq. (7) is equal to zero for r > if Sff(u) also has no poles in 
the lower half-plane. Since the spectral density satisfies the condition S(—u>') = >S(u/), 
its poles must be placed symmetrically with respect to the real axis uj = uj' . Therefore, 
the regression theorem can be fulfilled only in the case, when Sff(uj) is an analitycal 
function in the whole plane, i.e. it has no poles anywhere. 

Usually, to determine the properties of S/f(uj) under thermal equilibrium one uses 
the QFDT which relates the spectral density S xx (u) with the imaginary part of the 
susceptibility, Im{a x {u)}, by the simple relation: 

S xx {lo) = g(u)Im{a x (u))} (8) 

where the relating factor g(u) = hcoth(hu/2kT) is a universal function of frequency ui 
and temperature T only. Note, that the proof of the QFDT, given for the first time by 
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Callen and Welton [2] and in a somewhat other form in [3,4], is strictly valid on the 
real axis, uj = uj', only. In applying the QFDT, the expression (8) is usually extended 
to the whole frequency plane. The imaginary part of a x (uj) has a natural analytical 
continuation, defined as: Im{a x (uj)} = i[a{—uj) — a(u)]/2, however the relating factor 
g{uj) is continued into the complex plane as it is, without providing a sufficient physical 
justification. By assuming that such a continuation is valid and by comparing Eqs. (6) 
and (8), one obtains the main relation which is used to determine the properties of 
Sff(u>) under thermal equilibrium: 



where Im{A(uj)} = Im{a x (uj)} /[a x (uj)a x ( — uj)] is the imaginary part of the operator 
A t {} written in the spectral representation. This imaginary part is responsible for 
the relaxation processes of fluctuations in the system under test (i.e. it determines 
the relaxation law) and corresponds to the dynamical force of a friction written as 
fu = iIm{A(u>)}x u! . As follows from Eq. (9), the analytical properties of Sff(uj) 
are governed by two factors: (i) the poles of g(uj) coming from the QFDT, and (ii) 
the relaxation law of the system, i.e. the frequency dependence of Im{A(uj)}. For 
the regression theorem to hold, it is sufficient that the poles coming from g{uj) are 
compensated by the zeros of Im{A(uj)}. For example, this compensation takes place 
in the classical domain, where % — > 0, g{uj) — > 2kT/uj, i.e. there is a pole at uj = 0. 
To compensate this pole it is sufficient that Im{A(uj)} = uj(, where C is a constant 
independent of frequency. This case corresponds to the usual viscous friction used to 
describe the relaxation within the classical limit when Sffiyj) = 2kT( and the regression 
theorem is fulfilled. 

In the quantum domain, where h ^ 0, but the relaxation is still described by the 
same classical law, when fi = 0, Sff(uj) contains uncompensated poles at the Matsubara 
frequencies uj m = ifim, where fi = 2nkT/h and m = ±1,±2, ... This is the source of 
violation of the regression theorem, as demonstrated in [7] by calculating the right hand 
side of Eq. (7) for the damped oscillator. However, in any case this result can not be 
considered as a convincing proof of the QRT violation, since a relaxation process within 
the quantum description does not necessarily must satisfy the same laws as within the 
classical description. Thus, the validity of the QRT is directly related with the analytical 
behavior of the relating factor, g(uj), in the QFDT. Consideration of this problem is the 
task of the following section. 

The fluctuation- dissipation theorem. From the general principles of quantum statis- 
tics the correlation function of fluctuations of an observable x can be written in the 
operator representation as: 



while the corresponding function of the linear response is given by Kubo formula [3]: 



where ps is the density operator which describes some stationary state of the physical 
system under test characterized by Hamiltonian H$ with a set of eigenstates ip n corre- 
sponding to energy eigenvalues E n . The time dependence of both C xx (t) and a x (r) is 



Sff(u) = g(u)Im{A(u)} 



(9) 



C xx (t) = ^Tr{p s [x(r)x(0) + x(0)x(r)}} 



(10) 




(11) 
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determined by values of the operator x(t) taken at two different time moments t = and 
t = t only. Therefore, in the matrix representation it is: x mn (r) = x mn (0)exp(— iu mn r). 
To obtain from Eqs. (10) and (11) the spectral dependences of S xx (u>) and a x (u) appli- 
cable in the whole complex plane, we make use of the standard procedure of calculation 
of the Fourier transform in such a case: 



x mn (T)exp(iuJT)dT = 



i(uj - u mn - iA') i(u - u mn + iA") (u - u mn ) 2 + A 2 

(12) 

Here the integration interval is subdivided into two parts: (— oo,0) and (0, oo). Inside 
the first and second part each of the Fourier integrals converges separately to analytical 
functions of 10 in the lower {to" < 0) and upper (oj" > 0) half-planes, respectively, 
and diverges in the opposite half-planes. In the divergent regions the value of the 
integrals can be definded only as analytical continuation of the expressions obtained in 
the convergent regions. The direction of the shift of the divergent poles is determined 
by the sign of the iA-terms in Eq. (12). Under the usual assumption that C xx (r) = 
Cxx(-t), such poles of S xx {u) are placed symmetrically with respect to the real axis, 
i.e. A' = A" = A. When A — > the Fourier integral of Eq. (12) tends to a 5-function 
which we shall label as: 2ti5^{u — uJ mn ), where subindex A indicates that the analytical 
properties in the complex plane are determind in accordance with Eq. (12). Since 
a x( T ) = when r < 0, the integral in Eq. (12) is calculated in the interval (0, oo) only, 
and hence, his value corresponds to the second term in Eq. (12). 

By going in Eqs. (10) and (11) to the matrix representation and by applying the 
procedure of Fourier-integral calculations given by Eq. (12) one obtains: 



S xx (uj) 
hlm{a x (uj)} 



x 



|2 



mn 



Pn + Pm 
Pn Pm 



5^{u-u mn ) (13) 



where p n =< i^n\ps\^n > is the probability to find the system in the eigenstate ip n . 
As follows from Eq. (13), the analytical properties of both spectra S xx (uj) and a x (u) 
are determined by the same set of symmetrical poles at u = uj mn ± iA, each of which 
describes an allowed transition between the eigenstates of the physical system. The 
difference of each spectrum is that the weight represented by the terms (p n ± p m ) with 
which a pole enters into S xx (u) and Im{a x (u)} is different, i.e. both quantities differ 
in the value of the residue at the poles while the set of poles is the same. Therefore, 
to interrelate S xx (oj) and Im{a x (cu)} (what is the content of QFDT) it is sufficient to 
formulate a procedure to recalculate their residues in each pole. 

Since at thermal equilibrium p m / p n = exp(—huj mn /kT) is a universal function of 
the transition frequency iv mn only, the full set of all weight factors, which perform such 
a recalculation in the poles of Eq. (13), can be described by a function of the current 
frequency given by: 

1 -p[u) 

where it is assumed that p(cu)| w=Wm „ = p m /Pm i- e - one obtains the expression of the 
relating factor in Eq. (8). Such a recalculation of the residues, is always present in 
more or less explicit form in any derivation of the QFDT. Usually it is performed in the 
following way. The multipliers 1 ±p{uj mn ) are isolated from the double-sum sign in Eq. 
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(13) as functions of to, i.e. 



hlm{a x (u>)} 



= 2tt 



"1 +p(u) 
l-p(uj) 



^ \%mn\ 2 Pn&A{u — U mn ) 



(15) 



Such a transformation is usually justified either by assuming a quasi-continuous spec- 
trum of the physical system [2] , or by using the properties of the 5-function inside the 
double-sum sign in Eq. (13) when A — > [3,4]. Then the double sum multiplier common 
to S xx (oj) and Im{a x (uj)}, which describes the resonances of the system under test, is 
excluded from consideration. As a result, one obtains QFDT given by Eq. (8) where 
the spectrum of the resonances is not given in explicit form. 

At the level of Eqs. (10) to (15) there is no reasons of QRT violation, since the 
properties of fluctuations and relaxation are determined by the same set of resonances 
of the physical system. The last step in the derivation of QFDT is the transition from 
Eq. (15) to Eq. (8). Here a formal violation of the equivalence of such a transformation 
can take place when Im{a x (u)} has zeros at frequencies u> m , where the multiplier 
1 — p(uj m ) = 0. Going from Eq. (15) to Eq. (8), the same multiplier 1 — p(ui) appears 
in the denominator of g{oS) so that the r.h.s. of Eq. (8) will have in hidden form 
uncertainties of jj-type, which can be a source of additional "properties" absent in Eq. 
(15). Since in the original Eq. (15) S xx (u>) is a regular function at frequencies u> = uj m 
it must remain a regular function in Eq. (8) too. This means that all poles of the 
relating factor g(u>) in QFDT, and hence, the poles of Sff(u) (see Eq. (9)), are points 
with removable divergence, i.e. the divergence at the poles must be compensated by the 
zeros of Im{a x (u))}. Thus, the correct solution of the uncertainty present in Eqs. (8) 
and (9) means, in essence, the equivalence of QFDT and QRT from a physical point of 
view. 

Generalization of the fluctuation- dissipation theorem. Contrary to thermal equi- 
librium conditions, in the general case p m j p n is not a quantity suitable to perform a 
recalculation of residues in the poles of S xx (u) and Im{a x (u))}. Usually, physical sys- 
tems are characterized by some set of transitions at the same frequency iv mn = u>r, as 
it takes place, for example, in the case of an equidistant spectrum. In this case p m j ' p n 
loses its unique dependence on frequency and can take different values for each separate 
transition. Accordingly, in Eq. (13) we regroup the double summation over the states 
of the system in such a way to put together all possible transitions with the same res- 
onance frequency ujr. Then, we number the whole set of different resonant frequencies 
in the whole frequency range [— oo, oo], for example, by the subindex R = 0, ±1, ±2, ... 
so that ujr + i > u}r, where ujr = iv mn is the resonant frequency of the system for all the 
transitions which satisfy the condition u mn = uj m t n i at n ^ n' and m ^ m' . In so doing, 
it is asssumed that ujr = uj mn and uj-r = u nm correspond to two different resonances 
with frequencies related as uj-r = —ujr. As a consequence, the double summation over 
the states in Eq. (13) can be rewritten as: 



^2 \ X mn\ 2 {pn ± Pmj^A^mn ~ Uj) = 

m,n 



R 



^ ^ Pn\%mn\ 
_E m =En+KujR 



i ^ ^ Pm\%nm\ 

E n =Em—hujfi 



(16) 
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Here the prime sign over the sum symbol inside the square brackets means that, when 
summing over the states the only contribution comes from n and m states satisfying 
the conditions reported at the bottom of the sum symbol. The first and second prime- 
sum in Eq. (16) at ur > collect all the transitions accompanied, respectively, by an 
increase and a decrease of the system energy by Tilor. Accordingly, when ujr < they 
interchange their places. It is easy to see that by defining the function p(u) as: 

P( W )L= WH == Yl Pm\x nm \ 2 / Yl Pn\x m n\ 2 (17) 

E n =Em-hLo R E m =En+hui R 

Eq. (16) can be written in a form analogous to that of Eq. (15), where p(u) satisfies 
the condition: p(— u) = p _1 (u;) and is uniquely defined in the whole set of resonant 
frequencies, i.e. at u> = ujr. Thus, it is always possible to write a generalized QFDT in 
terms of Eq. (8) where the relating factor can be represented in the sufficiently general 
form given by Eq. (14) and p(u>) is given by Eq. (17). 

In conclusion, as a consequence of the nonequivalent transformation used to derive 
the QFDT [2-4], it is shown that the QFDT given by Eq. (8) contains an uncertainty of 
{j-type originated by the poles of the relating factor g(u) and the zeros of the imaginary 
part of the susceptibility at the so called Matsubara frequencies. A correct resolution of 
such an uncertainty means, from the physical point of view, that the spectral properties 
of both fluctuations and relaxation are driven by the same set of resonances of the 
physical system under test only. In the framework of such a formulation, the QFDT is 
equivalent to QRT and the former can be considered as a mathematical representation 
of the Onsager hypothesis in the frequency domain. Moreover, a relation analogous to 
the QFDT can be written for an arbitrary stationary state of the physical system, by 
supposing that its linear response satisfies the causality principle. 

In the framework of the Langevin approach, the QRT is fulfilled under the require- 
ment that the spectral density of the force Sff(uj), which drives the fluctuations of 
observables, is an analytic function in the whole complex plane. Since the causality 
principle simply claims that there must exist some relaxation process in the physical 
system, but it does not detail the relaxation law, this requirement can be considered as 
a kind of selection rule for the choice of the relaxation law of fluctuations in order to 
be in agreement with the QRT. It is evident that such a requirement will call for some 
revision of the Langevin approach in the quantum domain. 
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